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Abstract 

In order to study the properties of J/'ip (IS) in the deconfining medium, we extend our previous 
formahsm [Phys. Lett. B 656, 45 (2007)] on J/tfj suppression at mid-rapidity using the colour 
screening framework. Our formahsm is more general as the complete rapidity, transverse momen- 
tum and centrality dependence including J/iJj suppression at forward as well as mid-rapidity can 
be computed directly from it. Careful attention is paid to the role of the medium's proper time in 
determining the locus of the screening region where J/ip gets suppressed. Other important ingredi- 
ents in the calculation are bag model equation of state for QGP, the longitudinal expansion of the 
QGP fluid obeying Bjorken's boost invariant scaling law and non-sequential/sequential melting of 
Xc (IP) as well as ip (2S) higher resonances. Upon comparison with the recent data of PHENIX 
collaboration on J/ip suppression at forward and mid-rapidity regions, we find that our model shows 
a reasonable agreement with the data without incorporating any sequential decay mechanism of 
higher charmonia states. Furthermore, we observe a larger suppression at forward rapidity in our 
model which is again well supported by the PHENIX data and also gives a hint that a scenario 
based on directly produced J/tp^s is preferable. 
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I. INTRODUCTION 



The current experimental programmes at RHIC are focused towards the investigations 
of the properties of hadrons or hadronic resonances above the deconfinement transition. 
In particular, suppression of J/il) (or charmonium state) has been suggested as a poten- 
tial probe of the deconfined matter produced in the ultra-relativistic heavy-ion collisions. 
Lattice quantum chromodynamics (QCD) calculations reveal that the critical temperature 
for confined normal nuclear matter (HG) to deconfined matter (QGP) phase transition is 
Tc ~ 0.160 — 0.190 GeV if baryon chemical potential = [1| although there are con- 
troversies whether such a transition is a first order or only a crossover P, iS] . Matsui and 
Satz 0] first predicted that the binding potential of the cc pair into J/ ip mesons is screened 
in the presence of a QGP medium and J/tp states dissociate at temperatures for which 
the colour (Debye) screening radius of the medium falls below their corresponding cc bind- 
ing radius. A crucial assumption in favour of QGP is that charmonia once dissociated in 
the medium cannot recreated at the hadronization stage. Thermodynamics of equilibrated 
QGP does not allow cc to be created in abundance and hence there is a suppression in 
J/ip yields. The strength of the suppression depends on the binding energies of the quarko- 
nia and the temperature of the medium [4i]. Old lattice QCD simulations [1,0] suggested 
that the J/'ip (IS) could survive in QGP up to dissociation temperature Tjj ~ 2.1Tc and 
the higher resonances, such as Xc (IP) and "0 (2S) can still melt near Tc- However new 
QCD analysis [tI shows that the dissociation temperatures of all the relevant quarkonia viz 
J/ip, Xc, and ■0 are close to Tc. The J ftp and their higher resonances are usually created at 
the initial (pre-thermal) stage of heavy-ion collisions because of their large masses. Their 
small widths also make them almost insensitive to the final state interactions. Therefore, 
J/ip quarkonia can probe the evolution of the deconfined matter beginning from the early 
stage of collisions [8]. 

There are now high-statistics data available on the J/ip suppression obtained by NA50 
experiments 0, [l3, at SPS and by PHENIX experiments at RHIC [11]. One of the 



surprising result is the observation of a similar suppression pattern at both these energies 
involving a large difference in the energy densities produced between SPS and RHIC. The 
other important feature of the PHENIX J/ip suppression data from Au+Au collisions at 
the center-of-mass energy ^snn = 200 GeV at RHIC is that J/ip yield in the central 
Au+Au collisions is suppressed by a factor of nearly 4 at mid-rapidity and 5 at forward 
rapidity as compared to that observed in p+p collisions scaled by the average number of 
binary collisions. The absolute suppression of J/^p in heavy-ion collisions can arise due to 
initial state effects (e.g., cold nuclear matter effects) and/or final state effects ( e.g., colour 
screening arising due to possible formation of hot QGP). Two types of models have been 
proposed in order to explain the similar suppression pattern observed for J/ip at SPS and 
RHIC energies. In first type of models, one considers that Xc and ^ states evaporate shortly 
above Tc while J/tp will do so if the temperature rises above 2.1 Tc or 1.2 Tc according to 
old [1] and new [3] estimates, respectively. Here one should emphasize that the new estimate 
for To for J/^/J poses an additional problem because then any enhancement in J/ip yield 
at RHIC energy would become difficult to explain. In second type of models, it is assumed 
that J/ip yield will result from a balance between annihilation of J/ip due to thermal 
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gluons 1^, 14] along with colour screening , |l^ and enhancement due to coalescence of 
uncorrelated cc pairs 17, 18, [l^ which are produced at RHIC energy 20, 21]. However, 
recent PHENIX data do not show a fully confirmed indication of J/ip enhancement except 
for the fact that {pj) of the data and shape of rapidity-dependent nuclear modification 
factor RAA{y) [13] show some characteristics of coalescence production. 

Cold nuclear matter ( CNM) effects such as nuclear absorption, shadowing and anti- 
shadowing 22] are also expected to modify the J/tp yield. CNM effects due to the gluon 
shadowing and nuclear absorption of J/ij: at the RHIC energy were evaluated from the 



J/il) measurement in d+Au collisions at RHIC [23l . |24 ]. Many pioneering workers have 



analyzed the SPS data [Ul] on J production (with regard to possible QGP effect) by em- 



ploying ideas such as colour screening model [la, |25|, |26[, nuclear absorption model based 
on cut-off energy density etc. Chu and Matsui approach was originally designed to 
study the pT dependence by parameterizing the energy dependence in the transverse z = 
plane but they did not include the bag pressure term B consistently. Karsch and Petronzio 
adopted essentially the same philosophy as Chu and Matsui without using bag constant 
and treating length of the plasma cylinder as a free parameter. Ruuskanen and Satz gen- 
eralized this method also to include rapidity dependence by proposing a factorized formula 
PiPT, y) ~ Upt) giy) for the fraction of deconfined cc pairs and studied its effect schemati- 
cally. Blaizot and Ollitrault 27] extended his nuclear absorption formulation [2^ by assum- 
ing that no J/V' survives if produced in a region where the energy density (proportional to 
number of participants) exceeds a cut-off value. They focused only on the centrality depen- 
dence without considering hydrodynamical expansion of the plasma. The current PHENIX 
data on the centrality dependence of J/V' suppression have already been analyzed in several 
models such as comover model fipl. statistical coalescence model S^, 31], kinetic model 17], 
statistical hadronization model[33] and QCD based nuclear absorption model [s^]- None 
of these models gives a satisfactory description of the present experimental data. In par- 
ticular, we do not find any single mechanism which can be used to explain the complete 
rapidity dependence of the J/V suppression. Chaudhuri has attempted to explain the sup- 
pression at forward-rapidity [3j] in a QGP motivated threshold model supplemented with 
normal nuclear absorption. However, this appears more like a parameter fitting because 
the parameters used in the analysis (e.g., QGP formation time r = 0.06 — 0.08 fm/c) do 
not convey any physical meaning. Moreover, in order to accommodate rapidity dependence 
in the formulation, additional three parameters are used in the analysis e.g., J/ij} nuclear 
absorption cross-section, threshold density and its smearing factor A to explicitly depend on 
the rapidity variable. Similarly Gunji et al. [ssl ] have recently used a hydro+J/'i/' model in 
which QGP has been considered as undergoing (3-|-l)-dimensional hydrodynamic expansion 
but it has been used to explain the mid-rapidity data only without explicitly incorporating 
J /if) formation time. Recently, we analyzed [15,] the mid-rapidity PHENIX data of Jjip sup- 



371 ] ■ using 



pression, normalized by contribution due to cold nuclear matter effects [35|, |36 
a modified colour screening model of Chu and Matsui [3] by improving substantially the 
basic theme of refs. 16|, l26|]. More precisely, we had taken a bag model equation of state, 
(l-l-l)-dimensional hydrodynamic evolution for QGP as well as additional time dilatation 
effect for charmonium as other ingredients in our model. We found that the model described 
the centrality dependence of J/ip suppression at mid-rapidity quite well 15]. 



3 



In this paper, we generahze the above geometric formulation 15|] in order to incorporate 
the rapidity dependence by explicitly introducing rapidity variable in a consistent manner. 
Here again we consider (l+l)-dimensional hydrodynamic evolution of the QGP and the 
possible sequential dissociation scenario [SQ] for the charmonium excited states. Our theory 
emphasizes the role of the medium's proper time in finding the locus of the screening re- 
gion and also shows how this locus explicitly depends on the J/ij: rapidity. We predict the 
centrality (i.e., impact parameter or number of participant nucleons Npart) dependence of 
the J/V' suppression in Au+Au collisions at mid-rapidity as well as at forward rapidity and 
compare with the recent experimental data of PHENIX collaboration at RHIC (normalized 
again with respect to CNM contribution). We notice that our model reproduces the main 
features of the data quite well. We also find that the present model gives more J/^jJ sup- 
pression at forward rapidity than that at mid-rapidity and this fact is again in quite good 
agreement with the recent PHENIX data. Results are reported corresponding to old and 
new [3] set of values of the dissociation temperatures of the quarkonia. 



II. FORMULATION 



Although our basic theme is similar to that described in refs. 15, la, 38], yet for the sake 
of convenience, we briefly recapitulate it below and also point out the important differences 
at ap pro priate places. For a QGP with massless quarks and gluons, the bag model EOS 
gives 



Mm- 



e = aTyci + B ; P = aT^ - B a {T^ - T^); (1) 

where at every time-space point x = {t,x), e is the energy density, P the pressure, Tc the 
critical temperature, Cg the square of velocity of sound and B is the bag constant. The 
coefficient a is obtained by assuming an ideal gas of massless quarks and gluons with zero 
chemical potential and includes the number of degrees of freedom as a factor since the 
-7r^/90 comes from integrating the thermodynamic distributions over momentum and is pro- 
portional to a value of Riemann-zeta function. It should be added here that the bag model 
EOS exhibits a first order phase transition between the QGP and the HG phases which 
is not well supported by the recent lattice QCD simulations [5, Even though the bag 
model EOS has often been used for hydrodynamic calculations, it should be understood as 
a convenient way of parameterizing some features of the EOS with a rapid change of en- 
tropy density as a function of temperature in the transition region. For a QGP undergoing 
(l+l)-dimensional Bjorken's boost invariant expansion, the local thermodynamic observ- 
ables become function of the lateral coordinate r along with the proper timer = (t^-z^)^/^, 
so that the cooling laws become 

i = P = f^ = f-1 ; q^l + d (2) 

where the dimensionless symbols f = r/Ti^T = T{T,r)/T{Ti,r), e = (e(r,r) — B)/{e{Ti,r) — 
B) and P = {P{T,r) + B)/{P{Ti,r) + B), have been introduced for convenience with Tj 
being the proper time for initial thermalization of the fireball. 
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From (1) the pressure is seen to almost vanish at the transition point T^, i.e., in the 
hadronic sector. Hence on any transverse plane we choose the initial pressure profile 

P{Ti,r) = P{n,0)h{r) ■ h{r)^l^-^ e{RT-T), (3) 

where Rt denotes the radius of the cylinder. The power (3 depends on the energy deposition 
mechanism, and 9 is the unit step function. Clearly, our pressure is maximum at the center 
of the plasma but vanishes at the edge Rt where hadronization occurs. The factor P(Tj,0) 
is related to the mean pressure < P >i over the cross-section and to the corresponding 
average initial energy density < e >j via 

P{n, 0) = (1 + /?) < P >,= (1 + P){cl <e>,-qB}. (4) 



In the present work, we take the initial average energy density < e >j in terms of the 
number of participating nucleons Npart [36] (which in turn depends on the impact parameter 
b), given by the modified Bjorken formula: 



< e >i 



Axn \dyH 



; At = 7tR't, (5) 



where At is the transverse overlap area of the colliding nuclei, {(IEt / dyH)yu=Q is the trans- 
verse energy deposited per unit rapidity of output hadrons. Both depend on the number of 
participants Npart [40] and thus provide centrality dependent initial average energy density 
< e >j in the transverse plane. The term ^ is a phenomenological factor. The motivation 
and need of the ^ factor will be discussed later in Sec. 3 in conjunction with the self-screened 
parton cascade model. 

It is well known that cc bound state in a thermal medium feels a colour screened Yukawa 
potential and it melts at the dissociation temperature To (determined by quenched and full 
Lattice QCD simulations 0, Q]) which corresponds to the energy density and pressure Pg 
given by 

Td>T, ■ es = aTl)/cl + B ; Ps = aT^-B. (6) 

For any chosen fireball instant t and on the arbitrary z plane the contour of constant pressure 
Pg is obtained by combining the cooling laws (3) with the profile shape (4) to yield 

P . ^- + ^ = f-.. (7) 

P(n,Q)h{r)+B ^' 

Setting r = the maximum allowed tilde time t^q (during which pressure drops to Pg at 
the center) can be identified as 



Ps + B 

with P(T.j,0) read-off from (4). Thereby the said locus takes the more convenient form 



f ^ - B/jPg + B) 
f%-B/{Pg + By 



1-75tI =Hs{t)^^^ ^777^^- (9) 
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Our above result generalizes a similar expression derived by us 15| for the special case 



when the transverse plane passes through the origin and the proper time was t itself. 

Next, consider an interacting cc pair created at the early time ti ~ h/2mc ~ at the 
location {ri,(pi,zi) inside a cylinder of length Li. The precise value of Li is not known 
a priori since different physical arguments can give different results, among which a recipe 
exploiting rapidity will turn out to be logical as discussed in the numerical section. The cc 



pair has mass m, transverse momentum pT, transverse mass niT = ym"^ + Pti rapidity y, 
total energy = rriT cosh(y), longitudinal momentum pz = rriT sinh(?/), vectorial velocity 
^ = {pt + 'Pz)/p^ and dilation factor 7 = p'^/m. In the fireball frame the pair will convert 
itself into the physical J/ip resonance after the lapse of time t = jtf (with Tp being the 
intrinsic formation time) provided the temperature T < Tjy. At this instant the pair's 
transverse position r, its longitudinal position z and medium's proper time r are given by 

r = ri+VTt ; z = zi+Vzt (10) 
T={t^ -z^ fl'^e{t-\z\) ; t = -fTF. 

From the locus (9) we deduce the so called screening radius 

rs = Rt{1- hI/^{t)Y'H{1 - Hs{t)]. (11) 

In contrast to the earlier paper [3], the important role of proper time r and hence of the 
longitudinal velocity Vz must be noted in the definition of the screening radius r^, which 
marks the boundary of the circular region where the quarkonium formation is prohibited. 
Since r and hence Hs(t) decreases as \vz\ increases, it is clear that the radius of the 
screening region grows with growing rapidity. This fact will be utilized for interpreting our 
graphical results in Sec. 4. 

Due to the existence of the screening region the pair will escape and form quarkonium 
if \ fi + VTt \> rs implying 

cos(0i)>C ; C^M^lll^^, (12) 

2 ri r \vt\ 

where the role of the transverse velocity vt has become explicit. This trigonometric in- 
equality is equivalent to saying that —(j)max{ri,zi) < (pi < 4>max{fi, zi) where 



(l>max{ri,Zi] 



TT ifC<-l, 

vr - cos-i I C I if -1 < C < 0, 

cos'i \ C\ if < C < 1, 

if 01. 



In contrast to our earlier work 15| the symbol (pmax depends on both ri and zi. 



Finally, we must deal with the anomalous survival probability due to QGP effect namely, 
S{Npart,PT,y)- Suppose the overall probability distribution for the production of cc pair 
at general position with general momentum is factorized as P oc /(ri) g{pT,y), where the 
radial profile function is 




/(n)oc 1-75^ e{RT-n), (13) 
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with a. = 0.5 [16] and the momentum distribution g{pT,y) is left unspecified because it 
cancels out from the expression of the net survival probability in colour screening scenario 
defined by 



ymax 



which can be simplified as 



S = „2 r / ^^I'^MI" 7^ r / dzi(j)max{ri,Zl). (15) 

TtK^Li Jo I \rCT/ I J-Li/2 



In contrast to our previous paper [15( above expression contains a non-trivial integration 
on zi coordinate. It must be emphasized that this generic symbol S refers to a meson of 
given species and specified momentum. 

Often experimental measurement of S at given Npart or y is reported in terms of the pT 
integrated yield ratio (nuclear modification factor) over the range {pT)min < Pt ^ {PT)max 
whose theoretical expression would be 

jiPTUa. ^ gdir 

/ qdiry •'yPTlmin (^P,\ 

dpT ' 

where the superscript "dir" pertains to the directly produced J/ip mesons. 



In nucleus-nucleus collisions 41[, it is known [36t | that only about 60% of the observed 
J/ ip originate directly in hard collisions while 30% of them come from the decay of Xc and 
10% from the ip . Hence, the pT integrated "inclusive" survival probability (S™"^') of J/ip 
in the QGP becomes 

(5'"^') = 0.6 (5'^^'')^ + 0.3 + 0.1 {S"^'')^' . (17) 

The hierarchy of old dissociation temperatures [61 thus leads to sequential suppression pat- 
tern {4^, with an early suppression of V'' and Xc decay products and much later one for 
the direct J/ip production. However, with new dissociation temperatures 0] employing full 
lattice QCD (which almost mutually coincide) all the three species will show essentially the 
same suppression pattern i.e., the concept of sequential melting will not have any dramatic 
effect as will become evident later in Sec. 4. 

Now let us turn towards the numerical section of our work where both To sets will be 
employed. 



III. NUMERICAL WORK 



Table 1 gives the values of various parameters used in our theory and the following 
explanations are relevant in this context. The value Tc = 0.17 GeV is in accordance with 
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the lattice QCD results [i]. The choice cj = 1/3 is most common for free massless partons in 
an ideal gas, although for partons which carry thermal mass or interact among themselves, 
may be different like 1/5 (s^ . The selection (3 = 1 [l^ indicates that the energy deposited 
in the collision is proportional to the number of nucleon-nucleon collisions. The initial 
proper time Tj for QGP thermalization is taken as 0.6 fm/c in accord with our earlier 
work [isl] as well as [4^. Also, relevant properties of the various quarkonia in a thermal 
medium are displayed in Table 2. The pressure Ps corresponding to old and new dissociation 
temperatures [a, [3] for various quarkonium species are calculated from (6) and shown in 
Table 2. It is clear that, in going from J /'ip to Xc to ip , the Td values (and hence also 
Ps) decrease some what sharply in the old set [a| but rather slowly in the new set 0]. 
As regards the length Li of the primordial cylinder (in c = ?i = 1 units) the Lorentz 
contracted total radius of the Au+Au system gives a small value ~ 2Rau/100 ~ 0.14 fm. 
The alternative estimate based on primordial creation time is also equally small namely, 
Ti ~ l/2mc ~ 0.066 fm/c which, however, can be made more logical by remembering 
that the corresponding longitudinal locations of the created cc pairs are expected to be of 
the order \zi\ ~ ri sinh(|y|). Here y is the rapidity whose experimentally determined bin 
reported in PHENIX experiment is 1.2 < \y\ < 2.2. Hence the length Li of the cylinder is 
expected to lie in the range 0.10 < Li < 0.29 so that its approximate mid value Li 0.2 
fm can be used. 

Our numerical procedure proceeds through the following steps: 

(i) Before finding the centrality (or impact parameter) dependence of J/ip suppression it 
is necessary to know the initial average energy density < e >, in terms of the number of 
participants Npart- For this purpose, we extract the transverse overlap area At and the 
pseudo-rapidity distribution {dEx / dr]H)r}„=o reported in ref. [13] at various values of num- 
ber of participants Npart- These {dET/dr]H)riH=o numbers are then multiplied by a constant 
Jacobian 1.25 to yield the rapidity distribution {dEx / dyH)yu=o occurring in (5). 

(ii) As regards the choice of the ^ parameter in (5) two routes are open. In the first route per- 
taining to the old Td values following point is noted. Even though the Bjorken formula 
provides a qualitatively good estimate of the initial energy density yet, it under-estimates 
the same, which can cause the suppression of only Xc and -0 but not of J/ip. Hence, a 
scaling-up factor ^ = 5 has been introduced in (5) in order to obtain the desired < e >i= 45 



GeV/fm [4J] for most central collision. The relatively large values of our < e >i have 
the following justification : These are consistent with the predictions of the self-screened 
parton cascade model these also agree with the requirements of hydrodynamic simula- 



tion 



441 ] which fit the pseudo-rapidity distribution of charged particle multiplicity dNch/drj 
for various centralities already observed at RHIC, and these can cause melting of all the 
quarkonium species listed in Table 2. In the second route pertaining to the new set of Td 
values 0] the original Bjorken's formula for (e), (without any ^ parameter) is sufficient to 
cause the melting of all quarkonia species. The appropriate characterization of kinematic 
quantities in Au+Au collisions is presented in Table 3. 

(iii) Next, we calculate the time f^o for the pressure to drop to Ps at the origin and thereby 
deduce the screening radius rs with the help of (8,9,10,11) for J/^lJ mesons of given px and 

y- 

(iv) Next, the quantity C is computed from (12) which sets the condition for the quarko- 
nium to escape from the scrG6ning rGgion ancl the limiting Vcilucs of the (prnax 

{ri,zi) are 
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constructed using equation written just below (12). 

(v) Finally, the survival probability S, at specified y and fixed pT but varying Npart is eval- 
uated by Simpson quadrature using (15) from which px integrated {S'^^^) is also deduced. 

(vi) The same numerical process is used to calculate the (5'°'*'') for all other higher quarko- 
nia and finally, by including sequential melting of higher resonances, total (5'*"'^') has been 
calculated from (17). 

(vii) Since we have calculated the theoretical survival probability of J/ip due to the QGP 
effect, hence, for comparison with the actual experimental data at RHIC, following proce- 
dure is adopted. The experimental nuclear modification factor Raa [13] is divided by the 



contribution due to CNM effect m, \22\ 

so as to yield the experimental value of 

survival probability S^Qp due to QGP viz 

qexp _ Raa /-.qn 

where Raa is the standard nuclear modification factor and Raa^^ ^ contribution to Raa 
originating from CNM effects constrained by the data of d-\-Au collisions. It should be em- 
phasized that deuteron-nucleus dA data are interesting both to fundamentally understand 
the issues of quarkonia as well as CNM and also to separate these effects from hot nuclear 



matter. Often the CNM effects are calculated by using theoretical model calculation [23] 



that include either EKS [46|] or NDSG [47|] shadowing models for the parton distribution 
functions : in each case an additional suppression associated with a (Jbreakup of cc pairs 
is also included. These methods of calculating R'aa^^ thus model dependent. Uncer- 
tainties associated with these theoretical schemes can be avoided by using the so called 
data driven- method 37 1. In the data-driven method, the Raa^^ parameterized in terms 



of RdAi^y^b) as a function of impact parameter b and mutually opposite rapidities ±y 
summed over the number of nucleon-nucleon collisions Ncoii calculated by Glauber model. 
These three methods are employed in ref. |2J] and also Raa in Au-|-Au as well as Cu-|-Cu 
collisions are compared with corresponding i?^^^- 



We take Raa^^ from ref. 24 1 calculated by data-driven method. Finally by using 
(18) our input data for the survival probability are generated. Of course, for a confirmed 
inference regarding this issue, error bars on the said ratio must be known. For convenience 
the mathematical expression for the error bar in S is recapitulated in the Appendix 
using individual standard deviation of the numerator and denominator. Unfortunately, the 



PHENIX experiments [13, ]2M d° ^^^^ what the correlation coefficient r between Raa 
and Rj^^^ is. Assuming these quantities to be highly positively correlated i.e., r = 1, we 
find 



S 



A r, A pCJVA/ 

"ft A A ^AA 



Raa R^T' 



(19) 



Now we turn to physical interpretations of our results. 
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FIG. 1: (a,b) The variation of integrated survival probability (S) (in the range allowed by invari- 
ant pt spectrum of J/i/; measured by PHENIX experiment 12, 12^ ) versus number of participants 
Npart at mid-rapidity. The experimental data are shown by solid rhombus with error bars. In the 
top Figure, solid and dash dotted curves represent predictions of our present model without and 
with sequential melting using new set of T]j values 7]. In the bottom Figure, same are shown but 
by using old set of To values In both Figures = 0.6 fm/c is employed. 
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2: (a,b) Same as Fig. l(a,b) but at forward rapidity [12 



TABLE I: Various parameters used in the theory 
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0.60 
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<S> vs y at fixed central ity: new T ^ set 
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FIG. 3: (a,b) Graph showing the theoretical variation of pt integrated survival probability i.e., (S) 
with respect to rapidity y at fixed centralities. Top Figure corresponds to the results using new 
dissociation temperatures for quarkonia 0] with and without sequential melting. Bottom Figure 
corresponds to the same as in top but by using old dissociation temperatures [6|. In both Figures 
thin light curves and thick bold curves represent calculation with sequential melting and without 
sequential melting, respectively. 
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TABLE II: Masses, formation times, dissociation temperatures [65 and screening pressures Pg of 
J/iIj, Xc and ip' . 







Xc 


^' 


m (GeV) 


3.1 


3.5 


3.7 


TF (fm) 


0.89 


2.0 


1.5 


Td/T, (Old) 


2.1 


1.16 


1.12 


Td/T, (New) 


1.2 


< 1-0 


< 1.0 


Ps (Old) 


9.391 


0.507 


0.387 


Ps (New) 


0.638 


0.098 


0.098 



TABLE III: Kinematic characterization of Au+Au collisions at RHIC at ^/snJj — 200 GeV [12] with 
Ti = 0.6 fm/c in Bjorken's formula when old (with ^ = 5) and new sets of Td (without ^ factor) 
values d, 0] are used. 



^part 


Rt (fm) 


< e >i 
^ factor 


(GeV/fm3) (With 
in old Td sets) 


< e >^ (GeV/fm3) (With- 
out ^ factor in new Td sets) 


22.0 


3.45 


5.86 




1.17 


30.2 


3.61 


7.92 




1.58 


40.2 


3.79 


10.14 




2.03 


52.5 


3.96 


12.76 




2.55 


66.7 


4.16 


15.69 




3.14 


83.3 


4.37 


18.58 




3.72 


103.0 


4.61 


21.36 




4.27 


125.0 


4.85 


24.38 




4.88 


151.0 


5.12 


27.37 




5.47 


181.0 


5.38 


30.52 




6.10 


215.0 


5.64 


34.17 




6.83 


254.0 


5.97 


37.39 




7.48 


300.0 


6.31 


41.08 




8.22 


353.0 


6.68 


45.09 




9.02 



IV. RESULTS AND DISCUSSIONS 



In this section we shall present/discuss our numerical results remembering the fact that 
the cooling laws in (2) and hence the radius of screening regions in (11) depend on the 
proper time r. In turn when the J/V' meson's trajectory is also considered in accordance 
with (10) we have to distinguish between the cases of mid-rapidity (vz ~ 0) and forward 
rapidity {vz > 0) separately. 

Scenario at mid-rapidity 
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Figure l(a,b) show the variation of px integrated survival probabihty {S) with respect 
to number of participants Npart at mid-rapidity corresponding to the new and old To sets, 
respectively. The three curves on each Figures correspond to experimental data S'^Qp 
(cf.l8; solid rhombus), our model calculation without sequential melting {SjJ^) (cf.l6; solid 

line) and that with sequential melting (5*"'^^) (cf.l7; dash dotted line), respectively. It is 
obvious from both Figures that (S) decreases with increase in Npart both experimentally 
and theoretically because of the increase in energy density with respect to Npart inside the 
screening region. In Figure 1(a) both {SjJ^) and {S^^'^^) are in qualitatively good agreement 

with the experimental data although quantitatively (S'*"'^') is slightly underestimated. The 
tiny difference between S'^^J^ and 5'*"'^' shows that the sequential melting mechanism does 
not have any dramatic influence here due to a very small difference in the dissociation 
temperatures of the quarkonia species in the new estimate for To. Let us now pay attention 
to Figure 1(b) based on the old To set where a parameter ^ = 5 was needed to scale up 
the initial energy density. The agreement between S'fj'^ and Sqqp is excellent while that 

between 5*"'^' and Sq^p is definitely poor. The reason behind this observation is that 

since the dissociation temperatures of J/ij: is far above that of Xc and -0') the latter two 
species just break-up before they can possibly decay to produce secondary J ftp's. In other 
words, the directly produced J/V^'s are responsible for the observed survival pattern with 
sequential melting not playing any significant role. Furthermore, the difference between 
our theoretical predictions with and without sequential melting is substantial due to large 
difference between the dissociation temperatures of quarkonium species in the calculations 
with an old set of Td values. 

Scenario at forward rapidity 

Figures 2(a,b) are also similar to Figures l(a,b) except that they pertain to forward 
rapidity. The main difference between Figures l(a,b) and Figures 2 (a, b) arises from the 
fact that at fixed Npart but non-zero rapidity {vz > 0), (S) decreases with increasing rapidity. 
This is understandable because, as pointed out beneath (11) the radius of the screening 
region at higher rapidity is larger than that at lower rapidities. Therefore, for a J/tp meson 
of given px, the traversal time will also be correspondingly larger implying that the predicted 
survival probability must decrease with increasing \y\. Also, in Figure 2(a) both the curves 
of {SjJ^) and (5'*"'^') agree very well with Sqqp and reasonably well in Figure 2(b). Hence 
based on the analysis at forward rapidity, it is not possible to make a preferential choice 
between the {Sfj^) and (S*"^') models. 

Detailed rapidity dependence 

Figures 3(a,b) show {S'^") and (S*"'^') plotted versus rapidity y at fixed centralities. 
The chosen centralities are 0-5%, 20-25%, 40-45% and 60-65%. These centralities do not 
exactly correspond to the recent PHENIX measurement [12] on J /ip suppression. Rather 
they correspond to an earlier measurement [13] in which the relevant [dEx / dr]H)r^^=Q values 
versus Npart were given for insertion in Bjorken formula. It is clear from Figures 3(a,b) that 
{S) shows flat maximum around \y\ ~ and decreases for large y as explained above. 
This trend qualitatively agrees with the experimental data jl2i] but for doing quantitative 
comparison with the data, calculations will have to be done at exactly same centralities 
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reported in ref. [13] • Furthermore, the difference between the predicted values of {S) (both 
(S**^*^) and (S'*"'^')) between mid and forward rapidities diminishes in going from non-central 
to central collisions. This trend again agrees with the recent PHENIX data [12] • A satisfying 
feature of Figures 3(a,b) is that our model predictions for (S) are obtained by simply 



inserting the relevant values of the rapidity y in contrast to other model [3j| where different 



additional parameters are needed for describing the data in different rapidity bins. 

It should be noted that our model does not involve any free parameter although there 
are some parameters used in the calculation. We have assigned proper justifications to 
their values and these have also been used in other calculations e.g., in Chu and Matsui 
model fli| and in hydro dynamical models 



At this stage a comment is warranted on the uncertainties involved in the CNM ef- 



fects 2j| extracted by data-driven method and constrained by d-|-Au collision measure- 
ments. In view of the relatively large error bars in R^^^ reported in ref. 24], it is not 



possible to make any firm quantitative statements on any additional J/tp suppression in 
Au-|-Au collisions beyond that expected from cold nuclear matter effects (i.e., due to pos- 
sible QGP formation). We feel, therefore, that more firm conclusion about possible QGP 
formation can be drawn by analyzing results such as Figures(l-3) if more accurate measure- 
ment of RdA is done in future experiments. 



V. SUMMARY AND CONCLUSIONS 



Generalization of our previous work 15|, based on the J/V'+hydro framework involving 
(l-l-l)-dimensional expansion of Chu and Matsui [I^ model, has been done in the present 
paper in order to incorporate complete rapidity dependence of the J/ip suppression. The 
present work also takes into account the additional time dilatation effect for charmonium 
formation due to its motion along a general direction and this makes our model different 
from other calculations e.g., by Chaudhari 34] and by Gunji et al. [35*]. Our formula- 
tion for J/ip suppression shows explicit dependence on transverse momentum, centrality 
as well as on the rapidity and the predictions are consistent with the data. Our model 
has suitably incorporated the results from a recent lattice simulation on the dissociation 
temperatures 0, and formation time of J/ipi well as on the non-sequential/sequential 
melting of higher resonances. As expected, using new To values we find that within error 
bars, the predictions without sequential decays of higher charmonia states provide better 
results. We conclude from our results that sequential decay mechanisms usually invoked 
for explaining the similarity of J ftp suppression pattern at SPS and RHIC, is simply ruled 
out by our present model. The conclusion remains surprisingly valid even in the case of the 
old values of Tjj used in the calculation for different charmonia states. It appears that the 
presence of a large initial energy density in this case makes the survival of various charmo- 
nia states difficult one and consequently the feed down factor from higher charmonia states 
does not improve the situation. Thus our calculation over emphasizes the role of produc- 
tion and dissociation of direct J/tp in the QGP scenario. We have analyzed the centrality 
dependence of the J/^|J suppression data available from RHIC (normalized by contribution 
due to CNM effects) in terms of the survival probability versus number of participants at 
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mid-rapidity as well as at forward rapidity. Our results reproduce the main features of 
the PHENIX data since we observe more suppression at forward rapidity as compared to 
mid-rapidity region and this is consistent with the recent PHENIX data. It is creditable to 
notice that the same mechanism explains the suppression patterns observed in the entire 
rapidity regions. Further extension of this work, such as incorporation of (3-1-1 )-dimensional 
hydrodynamic expansion and the predictions of our model at LHC energies would provide 
additional support to our ideas. 

In conclusion, J/^ suppression can still be regarded as one of the most potential signa- 
tures for deconfinement. Recent PHENIX data on J /ij) suppression and theoretical Lattice 
estimates for the dissociation temperatures of J/ if) family have added more confusions re- 
garding the interpretation of the data. Does a low depletion of J/V' yields at RHIC support 
a much stronger direct J suppression (at dissociation temperatures close to Tc) and some 
other mechanism like abundant charm quark production thus creating more cc pairs ? More 
precise data for d-|-Au collisions are needed to draw a firm conclusion on this issue. We 
believe that the LHC measurements of J/if) yield in Pb-|-Pb collisions at 5.5 TeV will help in 
resolving some of these confusing issues like sequential screening scenario and dissociation 
behaviours of J /il) and/or T families. 
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APPENDIX : ERROR BAR CALCULATION IN S (cf.(18)) 



Let o be a random variate with mean (a), fluctuation 6a = a ^ (a), and standard deviation 
Aa = (5^)^/^. Similar quantities for another variate b are called (6), 6b and Ab, respectively. The 
covariance between a and b is denoted by Cab = {Sa 6b) = rAa A(, with the correlation coefficient 
r lying in the range — 1 < r < +1. The ratio S between the variates has the statistical properties 



As 
S 



whose special values read 



is 
s 



Sa 



b 



{C-i)Y/^ = {{^)' + {^)'-2^^r} 



1/2 



(Al) 
{A2) 



As 
S 



1/2 



{A3) 



according asr = +l, —1, 0, respectively. 

In the problem of J/ tp suppression the error bar formula (19) is obtained by identifying a = Raa, 
and assuming r = 1. 



b = R^^T' 
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